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Abstract 

We examine gravitational waves in an isolated axi-symmetric reflexion sym- 
metric NGT system. The structure of the vacuum field equations is analyzed 
and the exact solutions for the field variables in the metric tensor are found in 
the form of expansions in powers of a radial coordinate. We find that in the 
NGT axially symmetric case the mass of the system remains constant only if 
the system is static (as it necessarily is in the case of spherical symmetry) . If 
the system radiates, then the mass decreases monotonically and the energy 
flux associated with waves is positive. 
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I. INTRODUCTION 



The present work examines gravitational radiation in the Nonsymmetric Gravitational 
Theory (NGT) (for a recent detailed review see U). We probe the NGT asymptotic be- 
haviour in the wave zone using an exact solution. A complementary analysis based on a 
DeWitt style background field expansion has already been published, announcing the main 
result of this current work 0]. 

The motivation for this work is twofold. Firstly, the nature of gravitational radiation 
is an important physical question which must be addressed in any candidate gravitational 
theory. Secondly, the literature already contains several incorrect treatments of gravitational 
radiation in NGT 0-§f, and since more incorrect treatments are currently being published 
M, it is important that the record be set straight. The reasons why erroneous results were 
arrived at by these other authors are explained in 

In General Relativity (GR) gravitational radiation from bounded sources has been stud- 
ied not only through the linearized theory but also with the use of exact solutions. The latter 
was done for the general case of a bounded source in asymptotically flat spacetime [§] . It was 
found that confining the arguments to the axially symmetric case did not cause any essential 
loss of generality. Since even the relevant GR calculations are very tedious and the level 
of technical difficulty in the case of NGT increases considerably, we limit ourselves to the 
axi-symmetric case. The GR gravitational waves from isolated axially symmetric reflexion 
symmetric systems were studied in detail in |J . Since our treatment of the axi-symmetric 
NGT case is rather parallel, familiarity with this analysis is strongly recommended. 

Since NGT was introduced [l(J there have been few analytic solutions of the field equa- 
tions published. The exact solutions known to date include the spherically symmetric vac- 
uum case [|ll[] , the spherically symmetric interior case [|!2| , |l3f and Bianchi type I cosmological 



solutions with and without matter [14,T5|. This, at least in part, follows from the fact that 



deriving NGT field equations relevant for particular cases of interest is not as technically 
simple as may be suggested by its superficial similarity to the corresponding GR situations. 
Firstly, since the underlying geometry is non-Riemannian, neither the fundamental metric 
tensor g^ u nor the affine connection is symmetric. This does not constitute a serious prob- 
lem for the choice of the form of g^, since we can always assume that its nonsymmetric 
part takes on the isometries of the symmetric part, which in turn has a well defined GR 
limit. On the other hand, calculating the connection coefficients proves to be a tedious 
and time consuming exercise, independent of the method chosen. Secondly, the resultant 
formulae for the nonsymmetric connection coefficients are extremely complicated for all but 
the simplest forms of the metric, thus becoming unwieldy to use in the derivation of still 
more complicated field equations. 

The NGT quantities presented in this paper were derived with the use of symbolic alge- 
braic computation procedures. To this end, we have used the symbolic computation system 
Maple. 

In Section [TT], we briefly summarize the necessary fundamentals of NGT. Section [TTT| deals 
with the coordinate system and generalization of the GR metric to the NGT case. Then in 
Section [IV], we expand the metric in negative powers of a suitably chosen radial coordinate 
and analyze the field equations. The closing section contains our conclusions. 

Throughout this paper we use units in which G — c— 1. 
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II. NGT VACUUM FIELD EQUATIONS 



The NGT Lagrangian without sources takes the form: 

C = y/=gg^R llv (W), (2.1) 
with g the determinant of g^ u . The NGT Ricci tensor is defined as: 



1 

+ W % W ^ ( 2 - 2 ) 



and W^ v is an unconstrained nonsymmetric connection : 

W£, = Wfa + Wfa. (2.3) 

(Throughout this paper parentheses and square brackets enclosing indices stand for sym- 
metrization and antisymmetrization, respectively.) The contravariant nonsymmetric tensor 
g^ u is defined in terms of the equation: 

9^9.u = 9^9va = e (2.4) 
If we define the torsion vector as: 

W, = Wfa = \ (W; v - W^) , (2.5) 



then the connection: 



is torsion free: 



Defining now: 



= W£ + \siW v (2.6) 



r, = Ff m] = 0. (2.7) 



2 

r^r^ + r^r^, (2.8) 



we can write: 



R,u(W) = R, U (T) + ^W M , (2.9) 
where W[^ u ] = \{W^ V — W u ^). Finally, the NGT vacuum field equations can be expressed 



as: 



9^-9^-9^ = 0, (2.10) 



3 



W H ), = 0, (2.11) 



= I w m- ( 2 - 12 ) 



For the purpose of the analysis of Section [TV], it is convenient to decompose into standard 

symmetric and antisymmetric parts: R^), R^u], and then rewrite the field equation ( |2.12| ) 
in the following form: 

iV)(r) = o, (2.i3) 

%M( r ) = 0, (2.14) 
where we used equations ( |2.6|) , ( ^7|) and the notation: 

= R[H, P + fi H,,, + ^Wl,„- ( 2 - 15 ) 



III. THE METRIC 

Similarly to GR, the simplest NGT field due to a bounded source would be spherically 
symmetric. However, the NGT equivalent of Birkhoff's theorem (see e.g. [|I|]) shows that 
a spherically symmetric gravitational field in an empty space must be static. Hence, no 
gravitational radiation escapes into empty space from a pulsating spherically symmetric 
source. 

Following U we consider the next simplest case: the field which was initially static and 
spherically symmetric and eventually becomes such, but undergoes an intermediate non- 
spherical wave emitting period. Also, spacetime is assumed to be axially symmetric and 
reflexion-symmetric at all times. Because of the complexity of the field equations, we are 
forced to use the method of expansion to examine the problem. This approach, namely 
expanding in negative powers of a radial coordinate, was also used in the GR analysis || 
and naturally suits a wave problem. 

Due to the physical picture sketched above and to the fact that we are interested in 
the asymptotic behaviour of the field at null infinity, X, (in an arbitrary direction from our 
isolated source) polar coordinates x° = u, x = (r, 9, 0) are the natural choice. The "retarded 
time" u = t — r has the property that the hypersurfaces u = constant are light-like. Detailed 
discussion of the coordinate systems permissible for investigation of outgoing gravitational 
waves from isolated systems can be found in [|,[|. 

The covariant GR metric tensor corresponding to the situation described above is: 



( Vr^e 2 ? - U 2 r 2 e 2 ^ e 2/3 Ur 2 e 2 ^ \ 

e 2(3 

Ur 2 e 2i o -r 2 e 27 

-r 2 e" 27 sin 2 6 J 



(3.1) 
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with U, V, j3, 7 being functions of u, r and 9 was first given in |T6| . 



For any metric in polar coordinates, form conditions must be imposed in the neighbour- 
hood of the polar axis, sin# = 0, to ensure regularity. In the case under consideration 
we have that, as sin 9 — > 0, V,/3,U/ sin 9,j/ sin 2 9 each is a function of cos 9 regular at 
cosfl = ±1. 

In order to find the NGT generalization of the metric tensor (|3.1|) we require that the 
symmetric part of the NGT metric tensor be formally the same as the GR metric tensor. 
We then impose the spacetime symmetries of the symmetric metric onto the antisymmetric 
sector. This is achieved by enforcing £g g^u] — 0, where the Killing vector field, is 
obtained from £g . g{ IUJ ) = 0. The solution to this equation for the metric ( p.l|) yields the 
single Killing vector field £(i) = Cfyd^ = sin 2 Qd$. Imposing £^ = yields: 

£fi)%[H + gi^du^ + 0[3*]0 M £ ( 3 i) = • (3.2) 

This equation gives d^g^] = 0, but does not exclude any antisymmetric components. This 
is markedly different from the static spherically symmetric case where the above procedure 
excludes four of the six antisymmetric components. Without further simplification, the 
NGT calculation would involve ten independent functions and ten independent non-linear 
differential equations. This would constitute a huge increase in complexity from the system 
of four equations and functions found in the GR case. 

To make the problem tractable, we need to determine which antisymmetric functions can 
be set to zero. To accomplish this we note that the imposition of axi-symmetry splits the 
antisymmetric field equations ( 2.11 ), ( 2.14j) into two sets of three independent equations. 



(This can be seen directly from the block- diagonal form of the GR metric). The first set 
explicitly involves the three skew functions gf 01 ' , g' 02 ^ ,g^'- 

(V=gg [H ) =o (/* = (), 1,2), 



#[oi,2] = o . (3.3) 
These four equations are not independent due to the one identity: 

-ggW) =0 foi/ = 0,1, 2). (3.4) 







,U,fJ, 

The second set of four equations explicitly involves the three skew functions g^ 30 " 1 , g^ 31 ' , g^: 

-gg [iv] ) =o, 



iWi = • ( 3 - 5 ) 

These four equations are also not independent due to the one identity 

e^ vp R M , P = . (3.6) 

We may now choose to work with either set of three equations and three functions, noting 
that eliminating one set of three functions simultaneously eliminates the three corresponding 



5 



equations. We choose to work with the first set of functions and equations as they reproduce 
the usual static spherically symmetric solution when the u and 9 dependence is suppressed. 
The components that we discard correspond to the magnetic monopole-like solutions of 
NGT. These components can be discarded without loss of generality as an extension of the 
background field analysis of Ref. [Q] to include all six skew functions does not alter our 
conclusions [fnj. 

In view of the above, the NGT generalization of the metric tensor ( |3.1| ) is: 



9\u> 



\ 



-l e 2/3 _ f/ 2 r 2 e 2 7 




Ur 2 e 2 ^ + X 







e 2 ?-uj 





a 







Ur 2 e 2 ^ - A 


—o 


_ r 2 e 2 7 
















— r 2 e~ 27 sin 2 





(3.7) 



where u, A and a are functions of u, r and 9. The contravariant metric tensor is given by: 



,01 



.aw 



( -o 2 e~ 2 ^A g- 

20 „21 



V 



9 







9 







9 



02 
,12 



- e 



9 

^ - u 2 )Ae~ 




■27 







-2 P 2 7 



r- 2 e 2 ^ sin' 2 9 J 



(3.8) 



where 



r 2 sin 2 9 

A= , 

9 

g = r 2 sin 2 9\-r 2 {e A < 3 -(to- aU) 2 ) + ae 2 ^'(2X - a 
g to = [ r 2 (uj -aU- e 2fi ) + aAe" 2 ^, 



,02 



,12 



r' 



,2(3 



u)cre 



~ 2 ^A, 



e 2!i ) + e 2/3 ~ 2 ^(A - aVr- 1 ) + Xcue^A, 



[Ur 2 {a 

g^[(cu,a,X)^(-u,-a,-X)]. 



IV. THE FIELD EQUATIONS 

The analysis determining the form of the functions U, V, (3, 7, lu, A, a in our case is a 
natural extension of that given in detail by Bondi et. al. [f| for finding the forms of U, V, /3, 7. 
The requirement that the field contain only outgoing radiation at large distances from the 
source gives the form of 7: 

AM) , gM) 
7 = + — — + ••• 

Demanding that the solution have the correct static limit (or equilibrium configuration) 
leads to the following forms for U, V, (3 and 7 (unless otherwise stated, all coefficients in the 
general expansions are functions of both u and 9) : 
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V = r -2M + — + ..., (4.2) 

r 

B-\ Bo 

(3 = -± + -4 + . . • , (4.3) 



c C-±c 3 

T = ~ + ~ + • • • > ( 4 - 4 ) 

The skew functions, u;, A and a, are constrained by the requirement that the spacetime 
is asymptotically Lorentzian and admits inhomogeneous orthochronous Lorentz transfor- 
mations ||. This requirement demands that g^g^ — > as r — > oo. In our present 
coordinates this condition is satisfied if u;, A and a have the following forms: 

w = -^ + ^ + ..., (4.5) 

A = L + — + • • • , (4.6) 

r 

a = So + ^ + • • • • (4.7) 

r 

The functions M, c, C, Lq, and W2 are a H functions of integration. Bondi et. al. refer 
to c as the "news function" as it controls the form of the gravitational radiation in the 
symmetric sector. In an analogous way, L is the "news function" for the antisymmetric 
sector. Consistent with these identifications, we shall see that the solution reduces to the 
static, non-radiative case when both c and Lq are set to zero. The static limit tells us that 
only M and W2 can be non-zero when the system passes through its equilibrium position and 
these coefficients will be identified as the mass and NGT charge of the body, respectively. 

We begin our analysis of the field equations by considering the simplest set of field 
equations - the skew divergence equations (y/—gg^), v = 0. The 9 component of this set 
becomes: 

= S , U - 2{S ° C) ' U - Sl > u +..., (4.8) 
r 

and since a must equal zero when passing through equilibrium, Sq — S± — always. Inserting 
this information into the u component directly returns W\ = from the lowest order term. 
The remaining, r, component yields: 

= L cot 9 + L fi - W 2 , u + ■ ■ ■ ■ (4.9) 

At this stage of the calculation, it is not profitable to continue to work with the skew 
divergence equations, as the next orders also contain unknown coefficients from the symmet- 
ric functions. Somewhat surprisingly, however, we are already in a position to calculate the 
NGT charge of the body, and to prove that it is conserved. The NGT charge, I 2 , is defined 
by the Gaussian surface integral 



gg [0v] ) v Sx = - rw 2 sm9d9 =< W 2 >, (4.10) 
2 jo 




7 



where the brackets <> denote the angular average. The charge is conserved since 

l\ = ~ f W 2jU sin0d0=~ r(L sme), e d9 = 0. (4.11) 
2 Jo I Jo 

This follows from the fact that L must be regular on the polar axis. 

We now turn our attention to the set of field equations Ruu>) = 0. The affine connections 
that we require to construct these generalized Ricci tensor components are obtained by 
solving the system of 64 equations ( |2.10| ). The closed form expressions for the non-zero 
connection components are extremely lengthy. For brevity, we only display the connections 
expanded in inverse powers of r, and only to the order required for our present analysis. 
The list of expanded non-zero components is given in the appendix. 

Rather than provide an exhaustive list of the R^u) = equations, we shall simply exhibit 
the components to the order that we require to obtain the solution. We begin with the (rr) 
component which demands: 

this gives B\ = B 3 = and B2 = — c 2 /4. The (r9) component then gives 

Q = i± + u 1 + u l c + 2 + 2cc a te + (413) 

which yields U\ = and U 2 = — (eg + 2ccot#). Inserting these expressions into the (uu) 
component results in the important condition 

M )U = -c u 2 + - {cee + 3c e cot 6 - 2c) >u . (4. 14) 

We can now use the (ur) components of R{^, v ) = at next order to solve for V\, but first 
we must choose how we wish to write the function of integration contained in U3. Following 
Bondi et. al. we shall write U3 as 

U 3 = 2N + 3cc e + 4c 2 cot 9 , (4.15) 

where N is the additional function of integration (the reason that U3 is written in this 
way rather than as U3 = N, and that the second function of integration in 7 is written as 
C — |c 3 rather than as C, is that we can then identify iV and C as the dipole and quadrupole 
moments of the source, respectively). Now the (ur) equation at next order: 



can be solved to give 



= sin 2 9 (lie 2 - 5(c e f + 4c 2 c u - 6cc, 
+U3 sin 9 cos 9 — 8c 2 — 19cc # cos 9 sin 9 

+ sm 2 9(u 3 , e -2L 2 + 2V 1 ) , (4.16) 



Vi = -N t e - N cot 9 + c/ + 4cc e cot 9 

+ ic 2 (l + 8cot 2 #) + L 2 . (4.17) 
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It is interesting to note that V\ contains the first explicit difference between the symmetric 
functions found for GR and those found for NGT. Substituting the above results into the 
(99) and (u9) components of R^u) = produces the following auxiliary conditions on the 
multipole moments C, N: 

4C% = 2c 2 c u + 2cM + N cot 9 - N e - L\ , (4. 18) 

-3N V = Af a + 3cc u6 + 4cc u cot 9 + c u c e . (4.19) 

We see that the quadrupole moment of a source will decrease more rapidly in NGT than in 
GR due to the Lq term. 

For completeness, we shall now return to the antisymmetric sector where the skew di- 
vergence equation can now be used to obtain the additional relations: 

W 3 = 2S 2)d + S 2 cot9 , (4.20) 
S 2 , u = L x -2cL Q . (4.21) 

The only remaining antisymmetric field equation ([2.14|) : 

-R[01,2] = -R[01],2 + R[12},0 + R[20],l = 0, (4.22) 

gives to lowest order: 

(W 2td + L 1 + S 2tU ) u = 0. (4.23) 
This equation yields the additional relation: 

W 2fi = 2cU-2L x . (4.24) 



Analysis of the solution 

To demonstrate the physical interpretation of M, we consider the static limit. We can 
scale the function c for either one of the static periods to be c = (forsaking the 9 dependence 
of c limits us here to a static spherically symmetric system). We now remove the terms 
containing the functions iV and C since they correspond to multipole moments. Since there 
is no radiation during the static period, N = C = 0. The metric ( |3.7| ) tends now to its static 
spherically symmetric limit: 





= 1 - 


2M S 




- + 

r 






I 4 

''s 

2r 4 ' 


9(oi) 


= 1 + 
l 2 


9[oi] 


- s 




9(02) 


= 9[02] 


= 0, 


922 


= -r 2 


• 


933 


= -r 2 


sin 2 9, 



l_l_2MJi 

^•4 ty 5 



(4.25) 
(4.26) 

(4.27) 

(4.28) 
(4.29) 
(4.30) 



9 



where by M s , I 2 we denote the static limit of M and l 2 , respectively. 

Now a coordinate transformation from the retarded time u to the usual time coordinate 
t = u + r converts the above metric into the NGT static spherically symmetric metric [|I|,|lIf : 

, 2 l\ s . 2M S , , 2 2M ~ x , 2 

ds 2 =(l + p(l--^)dt 2 -(l--^) dr 2 

-r 2 (d9 2 + sm 2 6d(f) 2 ) , (4.31) 

and 



#[oi] 



Thus, the static spherically symmetric limit, M s , of the "mass aspect" M(u, 9) can only be 
interpreted as the mass of the system. Similarly, the static spherically symmetric limit, I 2 



of the "charge aspect" H^(w, 9) is identically the NGT charge of the system as shown by 
equations (|4.10|) and ( |4.11j ). 



If we define the mass m(u) of the system as the mean value of M(u, 9) over the sphere: 

- r M(u, 9) sin 9d9 =< M(u, 9) >, (4.32) 
2 Jo 



m(u) 



then c(u,9) completely determines the time evolution of the mass m(u). Integrating (|4.14|) 



and noticing that the second term does not contribute to the integral due to the condition 
that c be regular on the polar axis, we get 

dvn 1 r n 

m,u = ^r = -- c u 2 sm9d9. (4.33) 
du 2 jo 

Since we discussed here systems whose initial and final states are static, the physical inter- 
pretation of m(u) as the mass of the system is unambiguous. Analogously to the GR case 
the main result is as follows: 

The mass of an axially symmetric NGT system is constant only if the system remains 
static. If the system evolves in time (emits waves), the mass decreases monotonically. 

Since radiation is the only energy loss mechanism available to the system, the above 
proves that gravitational waves emitted by an axi-symmetric reflexion symmetric NGT 
source compatible with the metric ( |3.7| ) carry positive energy or, in other words, the flux of 
gravitational energy in NGT is positive. 



CONCLUSIONS 

We have proved that an NGT axi-symmetric system emitting gravitational waves has 
the usual GR-like asymptotic behaviour in the wave zone. The NGT contributions to the 
physical quantities decay rapidly with the distance from the source and the energy flux at 
spatial infinity is necessarily positive. 

While we concentrated on an axi-symmetric source to simplify computations, the validity 
of our result is not confined to this particular symmetry. This contention can be made 
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concrete by adapting the analysis made by Sachs ||, where it was shown that the axi- 
symmetric solution contains all the important features of any isolated, radiative system in 
GR. Sachs' result follows from considerations of the asymptotic nature of the spacetime, 
and these considerations are unchanged in NGT. Physically, we can argue that an axi- 
symmetric source provides a complete range of multipole moments to act as a source of 
gravitational waves, and thus provides a general test of the wave sector of any gravitational 
theory. Moreover, in the wave zone, the superposition principle may be used to construct 
the radiation pattern of any isolated body from a suitable sum of axi-symmetric solutions. 
Our result totally refutes the recent claims that NGT has bad wave behaviour and 
shows that aesthetically unappealing, phenomenological modifications to NGT |J are not 
necessary. 
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j- ■ 
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.0 
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Vi (Ui, e + c u ) - B hu - M, 
+ [f/ 2 (2(c - S^c,, + C/ M ) + [/ 2)(9 + c + c^x) + - 5 2jU - S 
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+ 
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(L Qfi + \uiL^r- 2 + 
AW 2 S 2 r~ 7 + 

r- 1 -cr- 2 + Q c 3 -3c)r- 4 + .., 

W 2 r' 3 + (W,/^ -W 3 + cW 2 - - S 2 ^ r~ 4 + 

U 1 + [U 1 (c-2B 1 ) + U 2 + c e ] r' 1 
U 3 + U 2 (c- 2Si) + 2C/i (Si 2 - S 2 - cSi 



r 2 3 / sin 2 0= E/i - cot + (C/ 2 + c - ct/i - 2B-JJ\ + 4c cot 0) r" 1 



+ 



C/ 3 - C/ 2 (c + 250 + 2U X (B x 2 -B 2 + cBi) 
-Ac(c fi + 2c cot 9)]r~ 2 + 



c ,« r 1 + ( ^c z c u - C iU | r 



+ ..., 



L (Ui - cot 6>)r 



-2 



(A20) 

(A21) 
(A22) 



(A23) 
(A24) 
(A25) 
(A26) 

(A27) 
(A28) 

(A29) 
(A30) 
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+ \W 2 (Ui cot 9 - 1 - c u ) + L^C/i - cot 9) 

+L (c e - 2B X XJ X + 2c cot 9 + cU x + U 2 )\ r~ 3 + (A31) 

rfis) = ^ + c ^ 2 + (3C - ^c 3 ) r- 4 + (A32) 

rfi 3] = - W 2 r~ 3 + [W 2 (2B 1 -c) + S 2 (U! - cot 9) - W 3 ] r~ 4 + (A33) 

r? 23) = cot 9 - cer- 1 + Qc 2 c „ - C>) r~ 3 + (A34) 

rf 23] = - Lor' 1 + [L (2B! - c) - L x \ r~ 2 + (A35) 
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and <?[23] = rF_i + F + . . ., where Di, E iy Fj are functions of u and 6. Repeating 
the calculation in ref. 2 with all six skew functions gives the following: At lowest or- 
der the skew divergence equation, {\f—gg^), v = 0, gives E 0:ll = =^ E = 0. The 
lowest order in i?{[02],3]} = gives F_i ;UU = =3- F_i = (as demanded by the quasi- 
periodic boundary conditions). The next order in the skew divergence equations then 
gives Ei jtt = =>■ Ei — 0. Putting all this into -R{[oi],3} = yields D 0iU = =3- D = 0. 
Since the two sets of three skew functions have remained uncoupled to this order, the 
solutions for gpi], #[02], <7[i2] are unchanged from those in ref. 2. At this stage our work 
is done since all leading terms (1/r terms in cartesian coordinates) have vanished. It is 
thus a trivial exercise to show that g[ 03 ], p[i 3 ], #[23] can play no role in determining the 
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